arXiv: 1509.069lOvl [quant-ph] 23 Sep 2015 


Solid state based analog of optomechanics 

Nicolas L. Naumann/ Leon Droenner/ Alexander Carmele/ 

Andreas Knorr/ Weng W. Chow,^ and Julia Kabuss^ 

^Institut fiir Theoretische Physik, Nichtlineare Optik und Quantenelektronik, 
Technische Universitdt Berlin, Hardenbergstr. 36, 10623 Berlin, Germany 
"^Sandia National Laboratories, Albuquerque, New Mexico 87185-1086, USA 

Abstract 

We investigate a semiconductor quantum dot as a microscopic analog of a basic optomechanical 
setup. We show, that optomechanical features can be reproduced by the solid-state platform, 
arising from parallels of the underlying interaction processes, which in the optomechanical case is 
the radiation pressure coupling and in the semiconductor case the electron-phonon coupling. In 
optomechanics, phonons are typically induced via confined photons, acting on a movable mirror, 
while in the semiconductor system the phonons are emitted by the laser-driven electronic system. 
There are analogous effects present for both systems, featuring bistabilities, optically induced 
phonon lasing and enhanced phonon loss. Nonetheless, the different statistical nature of the optical 
cavity and the electronic system also leads to qualitative differences. 
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I. INTRODUCTION 


Quantum optomechanics (OM) has been rapidly developing into a major research area 
over the last decadesi*^. Intense activities are fueled by interesting physics based on radia¬ 
tion pressure arising from momentum transfer between photons and optical surfaces. The 
phenomenon affects ultra-accurate measurements (e.g., gravitational wave detectioii^) and 
can lead to nanoscale motors and refrigerators^. 

This manuscript proposes and analyzes a different physical system for exploring a semicon¬ 
ductor (SC) analog of quantum optomechanics. It differs from the atomic-molecular-optical 
(AMO) approach to optomechanics in that radiation pressure is replaced by the electron- 
phonon interaction. The presence of electron-phonon interaction usually leads to heating 
in optoelectronic devices^. However, if optically addressing phonon assisted resonances or 
when properly modifying and harnessing electron-photon interaction in an acoustic cavity, 
the same combination of interactions can be useful, e.g., provide high-temperature single¬ 
photon generation^, optoelectronic self-cooling and memory in silicon photonics.There is in¬ 
teresting physics as well, e.g. enabling of strong-coupling physics, via a quantum-mechanical 
polariton-polaron superposition^*^. 

Several proposals of phonon cavities coupling an acoustic phonon mode to the electronic de¬ 
grees of freedom have been made^*^ and even though there are challenges such as anharmonic 
phonon decay or fluctuations in layer thickness, high-Q phonon cavities have been realized 
experimentally!^^—. Besides such a designed electron-acoustic phonon coupling, that enables 
the interaction with a selected acoustic phonon mode, many semiconductor nanostructures 
naturally inherit a strong effective single mode electron-phonon coupling mechanism, i.e. 
the coupling to longitudinal optical (LO) phonons!^*!^. At optical excitation, these systems 
may hnd applications as phonon lasers!^. 

A possibility of these new applications and physical effects comes from the semiconductor 
electron-phonon coupling mechanism being stronger than radiation pressure. Further, the 
solid state platform typically involves much larger phonon frequencies^*!^, especially in the 
case of LO phonons. This in particular enables a controlled selection of subsystems for effec¬ 
tive phonon lasing or damping only by the choice of optical input frequencies. The sketch in 
Fig. [IKa) depicts the basic cavity quantum optomechanics setup. Photons give momentum 
to the cavity mirrors upon reflection, which then changes the mirror position and affects 
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QD Phonon cavity 



Figure 1. (a) Basic setup for investigating cavity quantum optomechanics. (b) Proposed semi¬ 
conductor quantum optical approach to exploring acoustic phonon dynamics. The phonon cavity 
consists of materials of alternating acoustic impedance in analogy to an optical distributed Bragg 
reflectors (DBR) ^^d^ . 


ultra-high precision (atomic scale) measurements. The Hamiltonian has the form^ 


Hqm = hUca,vC^C + hflrab^b + iflEi 



( 1 ) 


Here, c (c^) is the cavity phonon annihilation (creation) operator and cUcav the cavity fre¬ 
quency. is the frequency of the mechanical oscillator and b (fet) the corresponding 
annihilation (creation) operator. The third term models the driving by an external laser 
and the last term is the coupling due to radiation pressure, where L is the length of the 
cavity, m the mass of the mirror. 

The sketch in Fig. [I](b) is our proposed alternate experimental configuration for the semi- 
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conductor system, with the Hamiltonian 
Hsc = + hucola^ 

+ ihEi 

+ ^Y1 + (^q “ ala^K + H-c. (2) 

q 

Here, or (a|^ or aj) are the annihilation (creation) operators of an electron in the valance 
band or conduction band, respectively. In this case, feq (6]^) is the annihilation (creation) of 
a phonon with the frequency cuq. In this case, the electronic transition is driven optically. 
Furthermore, the phonon modes will be treated as one effective phonon mode.— 

We want to emphasize some formal similarities of the Hamiltonians: In place of the photon 
number operator in the AMO case, there is the electron population operator. In conven¬ 
tional quantum optomechanics, the photons create or remove mirror motion. In our case 
[Eq.([2])], a laser-driven electronic transition creates or annihilates a phonon. The coupling 
constants are of course different. 

The Hamiltonians for both systems can be brought to the form of Eq. (|3]), where we in¬ 
cluded additional electronic transitions, which are excited optically, for additional quantum 
dots. We now introduce a common nomenclature for the elements of both systems for no- 
tational convenience and to emphasize the similar roles the components play for the effects 
studied below. Thus, in the optomechanical system the photon mode is of high frequency 
(HE), while in the semiconductor system a two level system, which models a quantum dot, 
constitutes the HF part. The low frequency (LF) component is the mechanical mode of the 
oscillating mirror in the optomechanical case, whereas in the semicondutor case the phonon 
mode inside an acoustic cavity, as discussed before, is the LF part. 

From the sketches [cf. Fig. [T], it may seem like the semiconductor version is more compli¬ 
cated than the AMO one. However, with the actual experimental setups, the exact opposite 
is true. Conventional quantum optomechanics can involve large and complicated setups, like 
the LIGO for detecting gravitational wave. But there are also microscopic optomechanical 
devices, which maximize the optomechanical coupling strength by small oscillator masses. 
To do this, a broad range of designs exists in addition to the classical Fabry-Perot-Setup, 
using e.g. wispering gallery modes or photonic crystals^. 

The semiconductor system is on a chip, and very much like a photonic integrated circuit 
with the possibility of incoherent pumping of the semiconductor active medium. From the 
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AMO side, connecting the two approaches increases enormously the relevance of quantum 
optomechanics, because of the widespread use of optoelectronic devices. For example basic 
quantum mechanical features, such as the formation of mechanical polarons in polariton 
optomechanics^ as optomechnical counterpart to the polaron formation in typical semicon¬ 
ductor cavity quantum electrodynamics can be expected in quantum hybrid systems. From 
the semiconductor side, making the connection to AMO quantum optomechanics gives guid¬ 
ance on what may be interesting or useful to explore, transferring predictions and techniques 
for optomechanics devices onto optoelectronic platforms. 

With this guidance, we developed an analogy between the optomechanical coupling and 
semiconductor quantum-dot system and choose as demonstration of similarity between the 
two approaches three exemplary effects, which can be found in both cases. The manuscript 
is structured as follows. In Sec. m we introduce the Hamiltonian, which can be brought to 
the same form for both systems. The hrst effect we consider in Sec. IIIII is the bistability 
as a prototypical nonlinear effect. Next, mechanical lasing, also known as phonon lasing, 
is discussed in Sec. CYl It is present in optomechanical^'^'^ as well as in the semiconduc¬ 
tor systemsl^. Finally, in Sec. |Vl the basis of laser cooling of mechanical motion will be 
discussed, which can be observed in optomechanical systems22. Theoretically, it can be 
understood as an enhancement of the effective damping of the mechanical oscillator—. We 
compare the realization of these effects in the above systems. 


II. COMPARISON 

We compare two different physical systems to explore mutual effects but also to elabo¬ 
rate on differences. On the one hand, we are interested in an optomechanical system, which 
consists of an optically pumped cavity mode, which is coupled by radiation pressure to a 
mechanically oscillating mirror. As a semiconductor counterpart of such an optomechanical 
setup we investigate optically pumped quantum dots (QDs), that are coupled to a single 
phonon mode. 

The second case may be realized in two ways: The phonon mode is either a single conhned 
acoustic phonon mode or an effective mode describing the collective motion of several optical 
phonon modes. 

Both systems can be described by very similar sets of equations of motion, which will be 
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shown in the following. Thus, even though the main difference is the statistics of the HF 
mode, some differences will only arise from the choice of parameters. 

Despite all of the above mentioned effects being present in both, the OM and SC system, 
their origin is of different nature and thus also their experimental realization. As we will see 
in the following, the crucial difference between the systems is the behavior of the component, 
that is coupled to the mechanical degree of freedom. In the OM, the mechanical mode is 
coupled to an optical mode, which is bosonic and can thus have an arbitrary number of 
excitations. In the SC system, the phonon mode is coupled to one or several QDs, which 
constitute an electronic and thus a fermionic system. 

In this manuscript the QDs are modeled as a number of two-level systems (TLS). Therefore, 
the SC system can only contain a limited number of excitations. But also in more complex 
electronic systems the number of excitations is in principle restricted, which is not the case 
for the optical cavity in the OM system. This will lead to some major differences between 
the systems. E.g. cycling processes, discussed in the context of lasing, can be affected 
more strongly by electronic lifetimes. However, we will show that within certain parameter 
regimes and a selection of quantities the SC system can approach the OM one for a large 
number of QDs. Then the possible number of excitations becomes very large, as in the OM 
case, which would strictly be the limit of an inhnite number of QDs. 

Hamiltonian — The basis of our analysis constitutes a Hamiltonian model. The Hamil¬ 
tonians can be written in the same form for both systemsi^i^^ 

H = nnh^b + hoop^p + hgp^p (h + + iKE . (3) 

The terms can be identihed with the different components of each system: The hrst term is 
the mechanical or the phonon mode in the optomechanical and the semiconductor system, 
respectively. They are modeled as a single harmonic oscillator with frequency D. Then b 
(fel) is the annihilation (creation) operator for one phonon (we will also call the mechanical 
excitations of the oscillating mirror phonons), so the commutation relation reads [b,b^ = 1 
in both cases. This constitutes the part of the system with low frequency (LF) in comparison 
to optical frequencies. The second term models the cavity mode or the QDs. As shorthand 
notation, we introduce the vector p with elements pp whose length depends on the system. 
In the OM system, there is only one element, namely the optical cavity, and thus pi = c. 
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The cavity is modeled as a harmonic oscillator oscillating at frequency oj. The annihilation 
(creation) operator for the cavity mode is c (cl) (it is also [c, cl] = 1). The commutation 
relations for the degrees of freedom of the QDs differ due to their fermionic nature. In this 
case, p is now a vector with Nqd elements Pj, where Nqd is the number of QDs coupled to 
the mechanical mode. In the following, we are interested in the case of Nqd identical QDs 
(inhomogeneous broadening is neglected). 

We can now identify pj = (pi = with the anti-commutation relations 

~ Here, the j ^^e the annihilation (creation) operator for an 

electron in the band p of the jth QD mentioned in the introduction, of the jth QD. Here, 
we see that either the optical cavity or the QDs constitute the high frequency (HF) part of 
the respective system, since they oscillate at optical frequencies. The third term gives the 
coupling of the two elements with the coupling strength g. The number of excitations in the 
system pip is coupled to the position of the LF mode q = {h + ¥)/\/2. The last term is the 
pumping term, which describes the external pumping of the coherence p by a classical held 
with frequency cul- Here, E is again a vector and consists of the individual pump strengths 
Ej of the QDs. 

Equations of motion — From the Hamiltonian ([3]) the equations of motion can be derived. 
For the expectation value of some operator O they are given by 

fid) = ^(p.d]) + 7 Tr( 6 £|p]p) + ^Tr(0£|i]p) (4) 

+ ^Tr(d£[(T^]p). 

C[X]p = 2XpX^ - PXp - pX^X, 

where damping and pure dephasing are taken into account in the Lindblad form. In the 
OM system the cavity mode is radiatively damped with 7 , which is caused by scattering 
to other modes, outcoupling and imperfections of the mirrors. The mirror is damped with 
K because of relaxation to the environment. In the SC system the damping rate of the 
excited electronic state is given by 27 and the damping rate of the phonon mode is k. These 
originate from couplings to not explicitly treated photon and phonon bath modes. The pure 
dephasing exists only for the SC system with optical phonons. 

From Eq. (j4]), a system of equations is derived. For simplicity, we assume coherent helds and 
factorize products of operators into products of expectations values. Finally, the excitation 
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of the QDs is assumed to be equal for all QDs and that they follow the same dynamics. Then, 
the dynamics of each QD is described by the same equations. In the OM case, there are two 
equations, while we need three for the description of the SC system. The additional equation 
occurs since, for the optical cavity system, the number of excitations can be factorized into 
the coherences: {&c) = {&){c). In contrast, (cr+j(T_j) cannot be factorized, as this would 
imply that one individual QD can be excited more than once. After transforming the HF 
component to a frame rotating at the frequency of the external pump laser cul, the equations 
of motion read (shorthand notations for the expectation values: P = B = (b), and, 

only for the SC system, U = {p\pi)) 

B = -{in + k)B -igNU 
P = (iA - 7 - 7 pd) P - ig{B + B*)P 
+ Pi {1-UtU) 

P = Pi(P + P*)- 27 P. (5c) 

Here, the negative sign in Eq. fl5bj) is valid for the SC system, while the positive one holds for 
the OM system. Then, pump term does not depend on the occupation on the HF component, 
and Eq. fl5cll becomes obsolete. Also, the occupation U = P*P. In the SC case N = NqD, 
while in the OM case N = 1, since only one cavity mode is considered. 

The equations for the coherences P and B are formally identical: The HF component rotates 
at the frequency given by the detuning from the external pump laser A = ui, — uj and is 
damped by 7 . The coherence of the HF component is coupled to the position of the LF 
component, which effectively leads to a shifted frequency. Additionally, it is pumped by the 
external laser. In the OM system, the pumping term is not limited, while the pump in the 
SC system is restricted, since the QDs are saturated at some point. Finally, the coherence is 
coupled to the excitation of the HF component. The last equation describes the excitation 
of the upper level to account for the fermionic behavior in the SC system and contains a 
term with pumping and a term with damping. 

As shorthand notation we dehne the expectation value for the number of excitations in the 
HF component as uhf = NU, which becomes for the OM system tthf = P*P- In the OM 
system this is the number of cavity photons, while it is the number of electronic excitations 
in the SC system. The expectation value for the number of excitations in the LF component. 


(6a) 

(5b) 


i.e. the phonon number, is dehned as nLF = B*B. The OM equations are also valid for the 
SC system, when an individual quantum dot is excited marginally, which was shown in the 
context of spin systems, modeled by TLSs^^. 


III. BISTABILITY 


The hrst example we will give to show the similarities of the systems is the bistable 
behavior. The optomechanical bistability has been considered.— The equation for the steady 
state can be derived from ([5]) by setting the derivatives to zero. Furthermore, the stability 
of the steady state can be investigated by linearizing the system of equations in vicinity 
of the steady state^^. The bistability in the phonon number is given in Fig. ^a.) for the 
OM system and in Figs. [2](b) and[2]^c) for the SC system with acoustic phonons and with 
optical phonons, respectively. In all cases a bistability can be observed. For understanding 
the origin of the bistability, we will have a closer look at the forces on the LF component, 
which are given by Ftot = (B - B*)/{y/2i). 

There are two competing forces acting on the LF component: The restoring force resulting 
from the harmonic potential Fh on the one hand and the force given by the coupling to the 
HF component F^. In the steady state, an equilibrium of forces Ftot = F^ + F^. = 0 must 
prevail. When referring to a steady state value rather than the variable, we will add the 
index s. In the steady state, the restoring force is the same in both cases and linear with 
respect to the displacement Qs = (Bg + B*)/ \/2 of the LF component 

Fh= Qs- ( 6 ) 

In contrast to this the force F^., originating from the coupling of the phonon mode to the 
electronic excitation, is a nonlinear function of the displacement of the LF component 


Fc,sc = 


Ng 


V2 


1 + 


y (^A-V^gqs) +(7+7pd)^ 
7+7PD 


( 7 ) 


The radiation pressure force in the OM case is obtained in the limit N = 1, ypo = 0 and 
Ef ^ [A - \/2gqsY + 7^ as 


-fc,OM — 


2Ef 


\/2 (A - \/2gqgf + 


( 8 ) 
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Figure 2. Stationary states of the phonon number ulf for (a) the OM system, (b) the SC system 
with acoustic phonons, and (c) for the SC system with optical phonons. Note the difference in 
pump strength for bistability to occur. The parameters are given in Tab. HI 


A bistability can only be observed, for detunings A < 0, because there lies the maximum of 
Fc, regardless of the sign of g. Then, up to three equilibrium positions may occur. 

The difference between the OM and the SC case originates from the fact, that the force 
acting on the phonon mode is affected by the number of excitations in the HF component. 
The radiation pressure force can be increased by increasing Fi, even if the OM coupling 
constant is small. For the SC system, the absolute value of Fc^sc is at maximum Thus, 
either the coupling constant g has to be large enough, or the number of QDs has to be 
increased. Otherwise, pumping the QDs more strongly does not increase the force acting 
on the phonons, so that it may not become large enough to exhibit a bistability. This 
difference can be traced back to the statistics of the HF component, since its excitation can 
be increased arbitrarily by pumping in the OM case, while it is limited too in the fermionic 
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Figure 3. The limiting behavior of the stationary states of the phonon number tthf for the SC 
system (a) with parameters for LA phonons and (b) with OM parameters. An increasing number 
of TLSs is shown: In the limit N ^ oo, the OM solution (dashed) is approached (parameters as 
in Etb)). 

case. 

In App. 1^ we show that the equation for the SC system approaches the OM one for N ^ oo 
and 7 pd = 0, when the pump strength is rescaled according to Ei ^A"qd. This is illustrated 
for the SC system with LA phonons in Fig. [3](a). When the number of QDs is much larger 
than the typical number of phonons, the OM solution is reached by the SC system. In this 
case the coupling constant g is large enough, so that already for A = 1 a bistability can be 
observed. 

The same transition is shown for OM parameters in Fig. [31(b). In this case there is only 
a bistability for a large number of TLS, since the coupling constant g is rather small. For 
both cases the ratio Aqo/u-lf ~ 10^ gives an estimate for the number of QDs, for which the 
SC system coincides with the OM one. Here, ulf is the stationary phonon number in the 
regime of the bistability. 
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With this, the major difference in Figs. |2] can be explained: In the OM case the phonon 
number is very high (in the order of 10®), while it is rather low in the SC systems, since 
a bistability with a small coupling constant can only be observed, if the radiation pressure 
force is increased sufficiently through a high intensity in the optical cavity. In the SC systems, 
the bistability occurs already for very few phonons, since the electron-phonon coupling is 
stronger. 

Note, that in the SC system with optical phonons also pure dephasing is taken into account. 
Then the system is in a regime, where for one QD no bistability can be observed, but ten 
QDs suffice to obtain a bistability [cf. Fig. [21(c)]. With the additional pure dephasing, the 
limiting behavior shown above is not valid. 

After considering the bistability, which is a shared stationary property of the systems and 
motivates the OM system as a limiting case of the SC system with a large number of QDs, 
we will consider dynamical effects in the following sections. 


IV. LASING 

Fig. |3] suggests that the differences between the systems are more quantitative in nature. 
In the following, we will consider phonon lasing and damping. These processes, will reveal 
more profound differences between the two setups. At first, we will examine the effect 
known as mechanical or phonon lasing: A high number of coherent phonons is created 
in the mechanical mode. It has been investigated for the as well as for the SC 

systen>i^. Here, we focus on the comparison of the two systems. 

Induced phonon creation, that is exploited for phonon lasing can be understood as an anti- 
Stokes process of at least second order. For a most efficient build up of phonon population, 
the pump laser has to be blue detuned with respect to the HF resonance by roughly a 
multiple of the phonon frequency nVL. Due to a higher order energy conserving process, an 
excitation of the LF mode is achieved, when exciting the HF component [cf. Fig. 0]. At this 
point a major difference between the OM and the SC system already becomes apparent: For 
the OM system the upper- and the lower level in Fig. IHmay represent two arbitrary photon 
number states with n — 1 and n photons. Thus, the process is not restricted by the cavity 
decay rate k. However, for the SC system, the excitation process can only occur, if the QD 
is decayed into its ground state. Therefore, the possibility of creating a phonon within this 
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Figure 4. Scheme of the mechanism leading to lasing. See text for explanation. 


anti-Stokes process is limited by the decay time of the QD excitation. 

In the semiclassical model phonon lasing manifests itself as a parametrical instability: The 
steady state of low phonon number [cf. Sec. IIII] becomes unstable and the system reaches 
a periodic orbit of high phonon number. This is shown in Fig. [5] for the OM system and 
the SC system with LA and LO phonons. For the OM and the SC system with LA phonons 
the amplitude of the oscillations is small in comparison to the number of phonons since 
the damping of the phonon mode is rather small [See tab. [Tl]. This way, new phonons 
are generated via the above described process, while phononic excitation is still present. 
However, in the SC system with LO phonons, the phonon lifetimes are typically very small, 
while the radiative lifetime of the QD excitations remains the same. The creation of LO- 
phonons, in particular coherent phonons is therefore not as efficient [Fig. 0] and the phononic 
excitation decays for the most part before a new phonon is created. 

In Figs. [6] we show the mean phonon number in the system as a function of the detuning 
and the pump strength for the different systems. In case of the SC system, the number of 
QDs is further illustrated to have a crucial impact on the phonon lasing properties, such 
as thresholds or maximally achieved phonon number. For the OM system, we observe a 
maximum of the mean phonon number at the phonon assisted resonance for lower pump 
strengths [cf. Figs. [6][a)]. Pumping more strongly involves higher order processes, leading to 
a maximum at higher detunings. For the SC system, we note first of all, that the maximum 
phonon number is roughly proportional to the number of QDs [cf. Figs. |6][b)-(d)], which is 
consistent with the fact, that the QD excitations have to decay before the generation of a 
phonon can take place, i.e. in order to realize a phonon laser cycle. The next aspect, that 
can be observed, is that the maxima for phonon lasing are at the polaron-shifted phonon 
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Lasing: Time evolution 



Figure 5. The time evolution for exemplary parameters [cf. Tab. |TI] is shown for the OM (a), the 
SC system with LA phonons (b), and for the SC system with LO phonons. 


assisted resonances for lower pump strengths: In the case of a single QD iV = 1, only the 
hrst phonon-assisted resonance involving the hrst order anti-Stokes process (as depicted in 
Fig. H]) contributes to phonon lasing. For a larger number N of QDs [cf. Figs. [6](c),(d)] 
also the resonance involving two phonons becomes visible. In relation to this, we observe a 
shift of the resonance-condition for phonon-lasing. This can be evaluated approximately by 
considering an effective system accounting for the second order process depicted above. This 
is shown in App. [B], which follows Ref.—. The detuning from the effective QD resonance can 
be approximated as Aes = —A — 2-^ — 'UnFf^- The resonance given by the detuning (hrst 
term), as dehned inlTTl is altered by self-quenching (second term) with the pump strength 
and by the polaron shift (third term) with the electron-phonon coupling strength. For more 
than one QD [cf. Figs. [6][c),(d)] the polaron shift is proportional to the number of QDs. 
However, the QD excitation saturates fast, so that the shift is approximately constant in 
the lasing regime. The characteristic shift with the pump strength originates in the case of 
the SC system from the self quenching. 
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Figure 6. The mean phonon number is plotted as a function of the pump strength and the detuning 
for the OM system (a), the SC (LA) system with 1 (b), 5 (c), and 18 (d) QDs, and the SC (LO) 
system with 10 (e) and 50 (f) QDs. 


2 

For the OM system, the effective detuning is Agff = —A — This reveals, that there is 

no self-quenching, but also a dispersive shift, which is proportional to the number of photons 
in the cavity. This is analog to the proportionality to uhf in the SC case. Due to the weak 
coupling strength for most OM systems, the dispersive shift becomes only important for 
very high pump strengths. 

Further, for larger numbers of QDs [cf. Fig. EJ^d)], the one and two phonon lines mix so that 
the maximum for phonon lasing not only involves a first, but also a second order phonon- 
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assisted process. Then the above approximation becomes invalid. 

The behavior of the SC system with LO phonons is similar to the one with LA phonons, but 
with lower phonon numbers due to the inefficient creation of phonons discussed earlier [cf. 
Fig. in^e),(f)]. Nonetheless, the phonon number is still approximately proportional to the 
number of QDs, so that higher phonon numbers can be achieved by increasing the number 
of QDs. The mean number of phonons excited by the second order process is in a similar 
order of magnitude as the number of phonons in the stationary state for direct excitation. 
As in the case of LA phonons, the polaron shift and the self-quenching can be observed as 
well. 


V. ENHANCEMENT OF DAMPING 

The last effect, we will consider is the enhancement of the damping of the LF component 
through the interaction with the HF component. This effect is closely related to the lasing, 
which we examined in the previous section. It is also the basic mechanism involved in op¬ 
tomechanical backaction cooling.— For a full description of cooling, thermal fluctuations 
have to be taken into account, which will not be done in this manuscript. Here, we employ 
a simplified analysis and set the focus on the enhancement of the damping, which is usually 
treated in terms of a linear response analysis of the system of Eq. (l5]).-‘2i 
The effective damping rate of the LF component, which gives the time the system needs 
to return to its stationary state, can be altered by the interaction with the HF component. 
At hrst we will investigate the increase of the damping in the OM and the SC system by 
evaluating the largest Lyapunov exponent for the LF part of the linearized equations of 
motion, as computed for the stability analysis in Sec. IIIIl — This enables us to determine 
the most persistent oscillations in the phononic subsystem. 

In Figs. [71 the largest Lyapunov exponent Kefr,L is shown for the different systems as a 
function of the detuning A and the pump strength Ei. Note, that for the cooling in the OM 
case different parameters than for the bistability in Sec. IIIIl are used. This is done, because 
the stationary states become unstable shortly after the bistability for the parameters in Sec. 
mil and the enhancement of the damping for the stable states is rather small in this case. To 
avoid this, a set of parameters with weaker coupling is used [cf. Tab. [IT]. For the SC system 
it is possible to use the same parameters as in the case of the bistability, since the damping 


16 



Enhancement of damping 


bO^ 

O 
a ^ 

a 

Pk 


1000 

750 


0 


OM 


xlO^ 


Deviation x 10 


-2 


(a) 

1 

12 

(b) 

1 


1 

9 ^ 


= 

X 


6 



- 


3 


1 

^ 


0 

1 

1 


-4 g 


bX) 

a; 


sc (LA) 1 QDs 


a 

a' 

(Ih 


N 

m 


90 
60 - 
30 - 
0 


I 


6 

5 

4- 

3 

H 2 

1 


EK 

0 ) 


bjO 

d 

q; 

Cl. 

a 


N 

m 

o 


bq 


S C (LA) 5 QD s g 

1= 


90 
60 - 
30 - 
0 


4- 

3 

Ul 


Etl 

CD 


SC (LQ) 10 QDs 



.3 -1.0 
Detuning A/D 



SC (LA) 18 QDs^ 


(f) 


'k 

hP 

3 
H 2 


1 

■ 

1 

1 .0^ 


2 ; 

1 

0.9 1 



E 

0.8 

0.7 -1 


-^Ui 

Deviation x 10^ 



Detuning A/D 


Figure 7. Here, the enhancement of the damping for the different systems as a function of pump 
strength and detuning is shown, (a) The damping in the OM system is enhanced by a factor of up 
to 1400 with the chosen parameters, (b) The analytical formula ([9]) is in good agreement in this 
case, since the deviation is only in the order of 10“^. (c,e,f) For increasing numbers of QDs, the 
width of the enhancement of the damping increases, but it saturates at a factor of about 6 . The 
analytical formula (d) predicts an enhancement with additional pump and with the number of QDs 
[cf. (jlOp ]. For the SC system with LO phonons, the damping cannot be enhanced significantly (g). 
Even though the analytical formula only predicts a very high enhancement (h). 
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Figure 8. Scheme of the mechanism leading to enhanced damping. See text for explanation, 
saturates before the relevant stationary state becomes unstable. 

For the OM case [cf. Fig. EJ^a)], we observe an enhancement of the damping with increasing 
pumping power up to a factor of 1400. In the SC case with LA phonons [cf. Fig. [71[c,e,f)], 
the damping is also enhanced but only by a factor of about six. This is to be expected, as the 
cooling utilizes the process complementary to lasing [cf. Fig. |H]; When the HF component 
is pumped, due to the detuning in the order of the negative phonon frequency, a phonon is 
absorbed. This process can again only be as fast as the decay of the QD excitation, while 
this restriction is not present in the OM case. In contrast to the lasing case, where more QDs 
could increase the maximal number of phonons, the damping cannot be increased by adding 
additional QDs. This is due to the fact, that the stationary phonon number is rather small 
in case of the SC systems (< 1). This way not enough phonons are present in the system, 
to make use of the mechanism depicted in Fig. [Hlwith multiple phonons simultaneously. 

In the SC case with LO phonons, the damping cannot be computed via the largest Lyapunov 
exponent, since the phonon damping is not the slowest decay in this case [cf. Tab. Hi]- Here, 
we evaluate the damping numerically by htting an exponential envelope to the time evolu¬ 
tion of the phonon position. The damping determined this way is called K,eS,N- Fig. [71[g) 
shows the damping in the SC system for 10 QDs. The enhancement is so small, since the 
decay of the QD excitation is the same as in the LA case, while the phonon damping itself 
is already very high. Thus, the relevant process can only influence the damping marginally. 

For small coupling strengths or pump rates, an analytical formula for the effective damp¬ 
ing can be derived, as the dynamics in vicinity of the steady state is that of an effective 
damped oscillator without any interactions: B^s = — {i^es + «^efr) B. Here, we also in¬ 
cluded an effective frequency Des, since the coupling in principle also changes the oscillation 


high 

frequency 
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frequency. For the parameters used in our cases, however, the effective frequency can be 
approximated by Note, that the effective treatment is only valid when the dynamics in 
vicinity of the steady state can be described by a damped harmonic oscillator. In general 
more complex dynamics are possible so that the analytical formula looses its validity. 

To derive the formula, an external force acting on the LF component is taken into account 
in the Eqs. (jS]). This force could be a stochastic force originating from thermal fluctuations, 
which is needed for the analysis of cooling. The resulting equations are linearized in vicinity 
of the equilibrium positions and the result is Fourier transformed. From the resulting set of 
linear equations a susceptibility as the linear response to the external force can be computed. 
This is then compared to the susceptibility of a harmonic oscillator to incorporate the effects 
of the HF component on the LF component in the form of the frequency and the damping 
rate of an effective damped oscillator. More details on the calculations and all definitions 
are given in the App. O For the OM system, the effective damping rate is 2 i 


«^eff,OM [Ul) — n — 




7 ^ + (cn + a) 


7 ^ + (u — A 


which is valid for C*C ^ 1. The formula for the SC case reads 


(9) 


«efr,SC 


= K 


2Ng‘^EiVt 
^2 _|_ 4^2 

Re (P,) A - Im (P,) (7 - 7 pd )1 ( 7 ' - A' - 


2Re {Pg) A 7 + Im (P^) (277 + 


, V / 4P2 

7r-f/l + yWl- ^2 + 472 


7r + (+ A j 


'^1+ [ujR- A 


2 P 7 X (10) 


and is valid for weak excitation. 

In both cases A = A + ^^ 2 ^^ 2 “ ■ This includes the dispersive shift in case of the stationary 
state. This reveals why the dispersive shift is only observed in the OM system: The dispersive 
shift is proportional to the stationary number of excitations in the system, which is small 
for the SC system, while it is significant in the OM case. Even though the formula for the 
SC case is more complicated some similarities can be observed. In case of weak excitation 
the formula can be simplified by approximating ur ~ cu, 7 ^ ~ 7 , and P ^ 0. The major 
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restriction of the analytical formula is, that it relies on harmonic dynamics. For both systems 
this approximation fails at some point, but since the coupling in the OM system is weaker, 
the region of validity is very large and the deviation is only small. This is shown in Fig. [D^b) 
for the OM and in Fig. [7|^d) for the SC system. The deviation of the analytical formula from 
the full computation is given: For the OM case, only a small differences in the order of a few 
percent are observed [cf. Fig. EJ^b)], while for the SC systems, the formula overestimates 
the effect of cooling drastically, especially for large pump powers [cf. Figs. [71[f,h)]. 


VI. CONCLUSION 

In conclusion, we studied similarities in the behavior of an optomechanical and a semi¬ 
conductor system with a single phonon mode, constituted by either a single acoustic phonon 
mode or an effective optical phonon mode. We showed that the mechanical mode of the 
mirror has the same function as the phonon mode. Also, the cavity resonance and the op¬ 
tical resonance of the quantum dot are analogous. The underlying model is a Hamiltonian 
formulation, from which equations of motion are derived. 

Even though the systems exhibit similar features, the underlying physics are different. In 
particular, as already discussed in the model section [TTl the high frequency component is 
bosonic in the optomechanical case, while it is fermionic in the semiconductor case. 

Both systems exhibit the nonlinear effect of bistability. Here, the optomechanical system can 
be understood as the limit of the semiconductor system with many quantum dots. This is 
the case, since the additional equation, needed for the description of the fermionic quantum 
dots, is not necessary for weak excitation of the quantum dots. 

Furthermore, the detuning of the system resonance from the pump laser can be used to con¬ 
trol the behavior of the phonon mode: In the blue-detuned regime, the number of phonons 
is strongly increased, leading to phonon lasing. In the red-detuned regime, the damping of 
the phonon mode is enhanced. 

In case of phonon lasing, the number of coherent phonons is very large for the optome¬ 
chanical system since the relevant processes are not limited strictly by the lifetime of the 
high frequency component. For the semiconductor system, this limitation is given, since 
the quantum dot can only be excited again after its excitation has decayed. The phonon 
number can be increased by additional quantum dots coupled to the phonon mode for the 
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semiconductor system. In case of optical phonons, however, their short lifetimes prevent a 
signihcant occupation of the phonon mode. 

The damping may be enhanced for the optomechanical system and the semiconductor sys¬ 
tem with acoustic phonons. For optical phonons an enhancement of the damping cannot be 
observed. This is attributed to the choice of parameters, especially to the short phonon life¬ 
time. For the optomechanical system, the damping may be increased strongly before other 
effects dominate and prevent an enhancement. For the semiconductor system, the damp¬ 
ing may only be enhanced weakly. Even additional quantum dots do not allow a further 
enhancement, which would be suggested from the observations for lasing. 
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Appendix A: Limiting behavior of the bistability 

In case of the bistability [cf. Sec. IIII] the necessity for an equilibrium of forces Fh + Ff. = ^ 
leads to a cubic equation with up to three real roots in the bistable regime. When using the 
relation , which follows from (l5a|) in the stationary limit, the cubic equation 

can be expressed in terms of the excitation of the HF component. Note, that Utot = P*P in 
the OM case and Dtot = NU in the SC case. For the OM system, the equation, which has 
to be solved reads 



(Al) 


The corresponding equation for the SC system reads 




2 



+ (7 + 7 pd )' Utot = NEl 


If we neglect pure dephasing and assume 1, i.e. a single QD is only excited weakly, 

the rule discussed in Sec. Illll for connecting the systems can be observed by comparing Eqs. 
M and dMD. 
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Appendix B: Effective system 


To gain further insight into the processes, which are relevant for lasing, an effective 
Hamiltonian is derivedl^. This takes only the second order process discussed in Sec. nYi into 
account [cf. Fig. 0]. 

This is done in two steps: At hrst, a unitary transformation is applied in order to transform 
the Hamiltonian into a rotating frame. Afterwards, another unitary transformation is used 
to eliminate interactions in hrst order. The rotating frame is introduced by 

//rf = UHU^ - ihUdtU\ (Bl) 

U = 

^ = hcULP^p. 

This leads to the Hamiltonian 

H-rf = — hAp^p + hgp^p (b + 1^^) + iKE (p^ — p), (B2) 

'-^ ^ ^ „ 

where, as above, A = ur — u. This way, the explicit time dependence of the Hamiltonian is 
eliminated. In the next step, the unitary transformation 

Hes = ^ Ho + Hi + [tS,Ho + Hi] + ^[tS,[iS, Hq]] (B3) 

is applied such that only interaction terms of second order are present. By choosing S = 
Q:p + /5plp6 + H.C., the hrst order terms can be eliminated, and the ehective Hamiltonian 
can be computed from 

H,ii = Ho + ^[tS,Hi\. (B4) 

This is achieved with o: = —and {3 = i^. The ehective Hamiltonian reads then 



Here, we shifted the energy scale appropriately. The hrst two terms are the energy of 
the LF component and of the HF component in the rotating frame. The third term can 
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be neglected in the OM case, where it constitutes a shift of the total energy scale of the 
Hamiltonian (VFom = E^). In the SC case, this term introduces the self-quenching and can 
be written as Hsc = 2 ^iPlPi- The fourth term describes the dispersive shift. It 

depends on the number of excitations in the system. The fifth term models two processes, 
which are energy conserving in the case of lasing [cf. Fig. 0]: The number of excitations in 
the HF component is increased (decreased) by the external field and a phonon is emitted 
(absorbed). The last term describes processes, which are energy conserving in case of cooling 
[cf. Fig. [8]: The number of excitations in the HF component is increased (decreased) by 
the external field and a phonon is absorbed (emitted). Depending on the scenario one of 
the two terms can be neglected. 


Appendix C: Effective damping analytical 

The formulas for the effective damping are derived in analogy to Ref.—, while we use 
the semiclassical equations as starting point. For the linearized system, this leads to the 
same susceptibility^. At hrst the system of equations ([S]) is linearized in vicinity of the 
steady state. In addition an external force F on the LF component is taken into account; 
B = —{i^l + K)B — iNgU + iF/'/2. Then, the real variables q = , p = , X = 

and Y = are introduced and the the deviation from the steady state is defined as 

= (l{P) ~ Qs- Here, the index s indicates the stationary value. When terms, which are 
nonlinear in the deviations, are neglected, a set of linear differential equations emerges. This 
can be Fourier transformed, leading to the sets of equations 


— iUjFq = VLFp — K,Fq (Cla) 

—iojFp = —VtFq — tiFp — QxXx ~ QyXy + Xf (^^b) 

—iujFx = XFy — 7-Fx -|- QY^q (Clc) 

—%G}Fy = — ~XFx — iXy — gxXq (Cld) 
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for the OM system and 


—iOJpq = VtPp — Kpq 

(C2a) 

—iuPp = —Vtpq — Kpp — y/2N gpjj + Pp 

(C2b) 

—iujpx = APy — xXx + QvPq ~ 2^/2EiPu 

(C2c) 

-iupy = -APx - iPy - gxXq 

(C2d) 

—iujpu = \/2EiPx — 2Pu 

(C2e) 


for the SC system. The T indicates the Fourier transform of the respective variable. The 
Fourier transforms are always functions of a). We introduced the definitions 7 = (7 + 7 pd), 
A = A + \plgq^, gx = V^gXs, and gy = V^gY^. 

This system of equations can be solved so that, the susceptibility as the linear response to 
the external force can be computed from Xq = x Xp- In the case of the OM system, the 
inverse of the susceptibility reads 


XoM (^) 


+ 2iuK — 2g‘^P*PsA 

o (y — iujY + A^ 


For the SC system, it is 


(C3) 


Xsc (^) = 


+ 2iujK — 2gEiN 


gxA + gy {iu - 7 ) 


n 


ioj — 27 (—iu) (1 — i?) + 7 + 2R'y) (7 — iui) + A^ 


(C4) 


Here, we introduced the dehnition R = ■ Differences to Ref.— for the OM system arise 

from differences in the definition of the parameters [cf. Tabs. [H [H] and the introduction 
of damping. In the limit D ^ k, which is valid for the systems we are interested in, both 
descriptions are equivalent. From the susceptibilities, the effective dampings ([9]) and flTOll 
can be computed, where we use 


V = 


1 

2 1 - (1 + ff 


7r = 7 


Ur = U 



y/l + R- 
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Table I. Parameter values for the bistabilities in Sec. EH 



Optomechanical^^ 

Semiconductor (LA)^ 

Semiconductor (LO) 

Detuning of pump and 

A = -2.6 X n 

GO 

1 

II 

<1 

A = -H 

LE component 

D = 27r X 10 MHz 

D = 556.6 GHz 

D = 55.3 THz 

Losses 

opt. cavity 7 = 27r x 14 MHz 

TLS 7 = 5 GHz 

7 = 5 GHz 


mech. resonator k = 27r x 50 Hz 

acoust. cav. k = 0.5 GHz 

K = 100 GHz 




7pd = 100 GHz^ 

coupling strength 

\g\ = 952.7 Hz 

\g\ = 197.5 GHz 

\g\ = 5.1 THz 

Number of QDs 

- 

1 

10 
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Table II. Parameter values for enhancing or decreasing the damping used throughout Secs. IIVI and IVl 



Optomechanical2i 

Semiconductor— 

Semiconductor (LO) 

pump rate 

El = 60 GHz 

El = 80 GHz 

El = 9.01 THz 


[cf. Fig. [5Da)] 

[cf. Fig. [5l[b)] 

[cf. Fig. [5][c)] 

Detuning of pump 

A « 

A « 

A ~ zbD 

(+: Lasing, —: enhancement of damping) 




LF component 

D = 27r X 10 MHz 

VL = 556.6 GHz 

D = 55.3 THz 

Losses 

opt. cavity 7 = 27r x 2 MHz 

TLS 7 = 5 GHz 

7 = 5 GHz 




7pd = 100 GHz 


mech. resonator = 27r x 50 Hz 

acoust. cav. k = 0.5 GHz 

«; = 50 GHz 

coupling strength 

g = 205 Hz 

\g\ = 197.5 GHz 

\g\ = 5.1 THz 

Number of QDs 

- 

1 [cf. Fig. EDb)], 5, 18 

10 [cf. Fig. [5][c)], 50 







